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ABSTRACT 


Rayleigh  Wave  Rejection  by  Optimum  Filtering  of  Vertical  Arrays 

By 

Wm.  C.  Dean 


Optimum  array  processes  such  as  the  maximum  likelihood 
filters  are  usually  derived  from  the  cross  correlation  matrix 
of  the  time  series  outputs  of  the  array.  With  a  vertical  array, 
however,  the  correlations  of  both  the  incoming  P-wave  signals 
and  Rayleigh  wave  noise  are  predictable  from  the  structure. 

This  paper  presents  the  theory  of  maximum  likelihood  filters 
in  vertical  arrays  which  can  provide  undistorted  estimates  of 
the  signal  or  the  various  Rayleigh  modes  with  the  other  modes 
cancelled  out. 

The  performance  of  these  optimum  processes  will  be  dependent 
upon  the  validity  of  the  underlying  assumptions  such  as  horizontal 
bedding,  known  instrument  calibration,  and  knowledge  of  both  P  and 
S  velocities  with  depth.  Examples  are  shown  of  the  performance  of 
the  optimum  filters  using  synthesized  data  where  the  underlying 
assumptions  are  satisfied.  Then  these  same  optimum  filters  are 
applied  to  data  from  the  UBO  vertical  array. 


I.  THEORY 

vertical  arrays  otfar  some  intriguing  advantages  over  other 
arrangements  of  seismometers.  Assuming  the  noise  to  be  composed 
of  Rayleigh  waves  and  perhaps  soma  mande-propagating 
and  assuming  the  signals  to  be  teleseismic  P-waves,  the  noise  and 
signals  will  be  recorded  by  a  vertical  array  with  some  rather 
special  characteristics. 

we  can  imagine  that  each  mode  of  Rayleigh  wave  noise  possesses 

a  random  character  common  to  noise  functions  in  general.  However 

for  esch  Rayleigh  wave  mode  the  response  versus  depth  and  frequency 

relative  to  its  response  on  the  surface  is  predictable.  Bcca 

the  relative  depth  variation  is  predictable,  the  vertical  array 

can  be  summed  to  cancel  the  Rayleigh  modes.  In  contrast  the  noise 

over  a  surface  array  is  more  unpredictable  and  more  variable  and, 

therefore,  more  difficult  to  cancel  by  array  summation. 

normally  we  compute  the  attenuation  of  each  Rayleigh  wave 

mode  relative  to  the  surface  continuously  with  the  depth  for 

1,2,3,  .t,5,o„  are  thus  calculating 

several  particular  frequencies. 

mode  attenuation  (and  sign)  versus  depth  with  frequency  as  a  para- 
metsr.  However,  we  could  convert  many  such  computations  into  at¬ 
tenuation  versus  frequency  with  depth  as  a  parameter,  naturally, 
the  particular  depths  we  choose  would  be  the  levels  where  a  seis¬ 
mometer  was  recording  in  the  vertical  array,  in  this  manner  we 
can  compute  a  frequency  response,  ,  where 


T-ABponse  at  the  roth  depth 
““  response  at  che  surface 


•or  the  hth  higher  mode.  These  Rayleigh  modes  must  obey  these 
frequency  relationships,  H^„(w) ,  not  only  in  an  average  sense 
integrated  over  time,  but  also  in  a  transient  sense  for  each  in¬ 
stant  of  time.  Otherwise  these  noise  modes  would  not  be  surface 

waves. 
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The  v/aveform  for  any  of  the  rnodes  at  depth  will  be  either  in 
phase  or  180°  out  of  phase  with  that  at  the  surface.  Hence  the 
filter  response,  H,  (u;)  ,  will  be  an  ever;.,  real  function  as  must 

IITII 

its  transform  in  the  time  domain.  In  other  words,  these  filters 
are  phaseless. 

In  a  similar  way  teleseismic  signals,  and  whatever  mantle 
P-wave  noise  may  be  present,  will  possess  very  definite  delay 
(i.e.,  filter)  relationships  bebveen  the  various  outputs  of  a 
vertical  array.  Let  us  assume  that  the  signal  waveform,  a(t), 
is  the  same  at  all  levels  and  further  that  the  surface  reflection 
coefficient  is  unity  and  the  reflected  waveform  is  the  same  as 
the  incident  waveform.  Then  if  the  signal  recording  on  the  surface 
is  2  s(t),  the  signal  recording  from  the  rath  level  seismometer  is 
given  by 

2  (t)  =  s(t  +  T^)  +  s{t  -  T  )  (2) 

m  mm 

where  T  is  the  one  way  travel  time  for  the  signal  from  the  rath 
m 

level  to  the  surface.  The  frequency  equivalent  of  the  above  t^^e 
equation  is 

Z  (ou)  =  Si'ii)  2  cos(a;T  ),  (3) 

m  m  ' 

Thus,  similar  to  that  applied  to  the  Rayleigh  modes,  the  filter¬ 
ing  applied  to  the  signals  is  also  phaseless. 

We  will  define  the  noise  mode,  n,  ,  as  follov/s: 

k 

n^^  =  the  signal  mode,  s  ,  plus  any  mantle  P-wave  noise 
with  velocity  and  reflection  properties  identical 
to  those  of  the  signals; 

=  the  fundamental  Rayleigh  mode 

n^  =  the  first  higher  Rayleigh  mode 

n^^  =  the  k-2  higher  Rayleigh  mode 


-2 


The  composition  of  each  of  the  vertical  array  outputs  is  given  by 


k“l 

;  S,  ( U;)  =  ^ 

ZjCt) 

• 

=:  T. 

k 

;  =  1 
m 

m 

«2k<"' 

• 

2  (t) 
m 

=  r 

k 

*"k<‘> 

;  Z  (uj)  =  ^ 
m  k 

\k<”> 

where  these  equations  can  be  ejcpressed  either  as  functions  of 
frequency, w  ,  or  time,  t.  For  time  equations,  the  operation  (*) 
implies  convolution  between  the  h's  and  the  n's.  The  first 
recording  level,  Zj^(t),  may  be,  but  is  not  necessarily  on  the 
surface  of  the  ground. 

In  matrix  form  we  have 


Z  =  HF.  (5) 

The  vector  of  vertical  array  outputs,  Z,  is  measured.  The  matri 
of  filter  relations,  H,  is  computable  from  the  structure.  The 
problem  is  to  determine  each  member  of  the  unknown  vector  of  noise 
modes,  n,  by  a  filtering  process  just  as  a  seismometer  would  record 
the  mode  on  the  surface  with  the  other  modes  absent. 

If  the  number, m,  of  recording  levels  is  less  than  the  number, k, 
of  noise  modes  present,  (m,k) ,  then  the  system  is  underdetermined. 

In  this  case  the  noise  modes  cannot  all  be  separated  from  each  other. 
If  m  =  k,  we  have  a  determined  system  with  just  enough  recording 
levels  to  find  all  modes.  If  the  number  of  recording  levels  is 
greater  than  the  number  of  modes  present-  the  system  is  overdeter- 
mined.  In  this  case  all  modes  ca.a  bo  separated  with  the  extra 
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measureinents  used  to  overcow^  the  effects  failure  of  ageumptions, 
errors  in  calibrations,  and  instability  of  the  optimum  filters. 


Each  of  the  m  input  channels  must  be  filtered  and  these  results 
added  to  produce  each  doaired  noise  .wmo  in  the  output.  Itius 
the  solution  leads  to  a  matrix  rather  than  a  vector  of  optimum 
filters. 

The  filtering  operation  ia  therefore 

GZ  =  GHN  *  IN  »  N  (6) 

where  I  is  the  identity  matrix.  The  solution  G  »  will  not 
do  since  we  especially  want  to  find  :.olution8  when  H  is  not  square. 
Consequently,  we  let 

T  -1  T 

G  =  (H  H)  H  (7) 

for  which  the  operation  GZ  gives 

GZ  «=  HN  *  IN  *  N  (8) 

as  desired. 

Our  ftiethod  will  b<2  to  solve  these  equations  for  the  optimum 
filter  matrix,  G(uj)  ,  in  frequency  for  several  reasons.  First,  the 
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well  log  analysis  yields  the  H(u))  filter  matrix  for  the  Rayleigh 
modes  in  terms  of  period  or  frequency. ^  Second,  not  all 
Rayleign  modes  are  present  at  all  frequencies.  For  example, 'nly 
the  fundamental  mode  can  be  expected  tc  exist  at  the  lowest  fre¬ 
quencies.  Finally,  the  solutions  of  the  matrix  equatioi\s  are 
independent  of  each  other  at  the  different  frequencies.  Thus 
the  inversions  of  the  matrices  which  are  not  expected  to  be 
larger  than  lu  x  10  are  readily  handled  on  even  the  smaller 
digital  computers. 
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II,  SYNTHETIC  EXAMPLE 


A  vertical  array  of  five  vertical  cotppcnent  seisrnometers 
existed  in  a  well  at  the  UJuta  Basin  Seisraologicnl  Observatory 
(UBSO)  .  Figure  1  shov’G  a  schematic  diagram  ef  the  seismometer 
placement  in  the  vertical  array  and  the  log  of  the  surrounding 
lithology.  The  Rayleigh  dispersion  curves,  also  sliown  on 
Figure  1,  indicate  that  the  fundamental,  first,  and  second 
higher  Rayleigh  modes  will  all  be  supported  over  the  signal 
frequency  range  of  0.5  to  2.0  cps.  Since  we  want  more  record¬ 
ing  levels  than  modes,  th  se  tnree  Rayleigh  modes  plus  a  signal 
mode  will  be  the  only  ones  crnsidered  to  exist  in  our  synthetic 
modeling  of  this  well. 

The  Rayleigh  wave  dispersion  analysis  of  this  well  log 
leads  to  the  mode-depth  frequency  responses,  H  (uj)  ,  defined 

Ivin 

by  equation  (1)  . 

Figure  2  shows  plots  ot  these  filter  responses  for  the 
first  three  Rayleigh  modes  relative  to  the  responses  at  the 
surface  of  the  ground. 

The  optimum  filter,  G,  solutions  nre  the  least  squares 
inverse  to  the  Rayleigh  signal  filter  matrix  H.  Thus  the 
matrix  products  of  G  and  H  gives  the  identity  matrix  as  shown 
in  equations  (6)  and  (8)  GH  =  I.  In  our  test  model  we  have 
assumed  the  signal  exists  with  the  same  size  and  amplitude  on 
all  traces.  In  other  words,  we  are  neglecting  ghosting.  We 
have  derived  optimum  filter  solutions  only  over  the  signal  range 
of  0.5  to  2.0  cps.  For  the  frequencies  below  0.5,  v;e  have  smoothed 
the  frequency  responses  to  zero  according  to  a  sinusoidal  gain 
function  for  both  the  H  and  G  frequency  responses.  We  have  smoothed 
these  responses  to  zero  in  a  similar  way  over  the  range  from  2  to 
5  cps. 
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in  our  computations  the  zeros  in  the  identify  matrrx  are 

,e.e  than  1  x  lo'®.  Below  0.5  and  shove  2.0  ops  the  OH  produc 

,tiU  gives  zeros  in  the  oft  diagonal  elements.  The  ma 

Issumed  hclow  0.5  and  shove  2.0  ops  is  (sH)  where  the  scalar.  . 

is  the  cosine  function  of  frequency.  In  the  same  way 

tion  assumed  for  the  optimum  rilters  in  this  range  -  • 

consequently,  the  matrix  product  of  these  frequencies  rs  gr 

2  (9) 

aG  aH  =  “a  I 

consequently,  the  orthogonality  of  our  optimum  filter  solutions  is 

o... 

In  order  to  de..ranstrate  tnat  our  uy 
indeed  correct,  we  created  some  synthetic  data  which  does  conform 

L  Che  well  log  analysis  we  have  assumed  for 
chose  four  independent  seismograms  from  ou*  digi  a 

filtered  each  so  thMt  the  main  energy  co 
„e  handpas  liter 

between  0.5  and  z.v  cps.  filtered 

,  on  the  surface  of  the  ground  and  filterea 

P.R  the  fundamental  mode  on  tne  surre 

•  y  the  H  matrix  to  produce  the  equivalent  fundamental 
according  to  the  H  matrix  ^ 

t  the  various  instrument  depths  in  the  UBO  well. 

^ae  at  the  ,,,ner  mode  and  similarly 

second  trace  was  identi^iea 

filtered  according  to  its  depth  behavior  determined  by  th 

.  •  The  third  was  identified  as  the  second  higher  mode,  and 
„,atrix.  The  thi  artificial 

ThG^e  Rayleigh  mode  traces  we  t- 

a  e  a  mixture  of  signal  and  Rayleigh  modes  expected 
signal  to  produce  a  mixtur 

from  this  well. 

Fiaure  3  shows  the  data  from  tne  tiv 

V  if  the  signal  only  were  present.  The  companion  traces  ^ 
array  if  the  g  (estimate)  for  the  signsx. 

Figure  3  show  the  optimum  filter  output  ^ 

q.  1  fhe  first  higher,  and  second  higher  modes, 
the  fundamental,  tne  first  g 
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immediately  below  each  optimum  filter  estimate  is  the  tr 
„  n»ae  that  a  aeis«,.etet  wouia  see  on 

othot^noaaa  .«e  not  pteaent.  The  aietottion  of  - 

aught  ana  can  he  attrlhutea  to  the  h.napaes  propert.es  ot  our 
optia.ua.  fitter  aoiutrona.  The  payieigh_^ae  outputs  shouia 
-n  fact  are  less  than  7  x  10  everywhere, 
nglre  4  is  Bia.iiar  to  Figure  3  except  that  the  aata  contains 
the  signal  plus  the  funaaa,ental  eoae.  Ms  can  he  seen,  the  optxa«™ 
Ultlr  outputs  agree  in  site  ana  «avefor.  with  the  true  surface 
rasponses  of  the  signal  ana  the  £unaaa«ntal  Eaylergh  a.o  . 

Figure  5  is  srnilar  to  Figures  3  ana  4  exc  ? 
contains  the  signal  pic  the  funaamental  ana  first  higher  ay  i 
aes  Again  all  three  ™oaes  present  are  reproaucea  in  size  ana 

wave  shape.  The  output  trying  to  ptoo 

which  is  absent,  has  no  value  greater  than  1  x  10  . 

Figure  6  is  similar  to  the  previous  three  excep  . 
secona  higher  moae  has  been  aaaea  into  the  vertical  array  aata 
;.gain,  the  outputs  for  the  signal  ana  the  three  Rayleigh  moaes 
are  reproaucea  in  size  ana  waveform  within  a  few 
expectea  results.  The  signal  trace  shows  better  than  a 
improvement  over  the  best  signal/noise  ratio  available 

vertical  array- 
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Fifjure  6.  Optimum  Filters  Estimating  Four  Modes 
(The  Signal  and  Three  Rayleigh  Modes) 
from  Synthetic  Data  with  all  Four 
Modes  Present. 


lit.  EXAMPLE  WITH  UBO  NOISE  AND  A  SYNTHETIC  SIGNAL 

We  have  seen  that  the  optimum  filters  perform  as  expected 
when  the  data  properties  match  the  dispersion  analysis  of  the 
well  log,  Wa  must  now  apply  these  solutions  to  data  recorded  at 
the  vertical  array.  To  concentrate  on  Rayleigh  m'^de  suppression 
and  reproduction,  we  have  avoided  complications  which  might  be 
introduced  by  deghosting.  Any  signals  recorded  by  \.his  array 
will  have  ghost  reflections.  Consequently,  we  must  add  a  syn¬ 
thetic  signal,  which  satisfies  our  identicd  size  and  waveshape 
assumption,  to  tiie  five-level  noise  output  of  the  array.  The 
first  five  traces  on  Figure  7  show  the  noise  background  and 
synthetic  signal  used  in  this  test. 

There  are  several  sets  of  signal  and  noise  mode  estimates 
shown  on  Figure  7.  All  of  the  optimum  filters  used  require 
inputs  from  all  five  levels.  However,  we  can  ask  for  estimates 
of  only  a  single  mode  or  for  estimates  of  several  modes  in  the 
filter  outputs.  We  designate  a  five  level  input,  two  mode  out¬ 
put  set  of  filters  as  a  5  x  2  set.  The  more  modes  a  set  of 
filters  estimate,  the  larger  the  filter  gains  will  be.  It 
follows  that  any  errors  resulting  from  incorrect  dispersion 
analysis,  instrument  calibrations,  or  non-parallel  layering  will 
be  amplified  much  more  in  these  filters  with  fewer  degrees  of 
freedom. 

The  outputs  on  Figure  7  are  from  a5xl,  5x2,  5x3, 
and  a  5  X  4  set  of  filters.  The  5x1  set  estimates  the  signal 
merely  by  averaging  the  inputs  from  all  five  levels.  These  out¬ 
puts  show  progressively  increasing  low  frequency  errors  on  the 
estimates  of  all  modes  as  the  degrees  of  freedom  decrease  (i.e., 
as  the  number  of  outputs  increase) . 
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Figure  7.  Several  sets  of.  optiTOuni  filters  are 
applied  to  UEO  noise  data  with  a 
synthetic  signal.  The  more  nradeo  a 
set  of  optinsum  filters  estimates,  the 
larger  the  filter  gains,  and  the  more 
errors  are  amplified.  These  data  con¬ 
tain  considerable  energy  below  0.5  cps 
for  which  the  filters  are  not  designed 
As  a  result  low  fregriency  errors  are 
large. 


V.  filter  solutions  do  not  agree  with 

we  know  thBt  the  optxmtn.  f*lte 

.  c  4"  0 s  bdov?  0«5  cps- 

the  well  analysis  for  f  .=4 iters 

tne  W«.>.  J  nather  the  optitnum  filters 

nnt  analvsied  below  this  frequency.  Fa 

nterpolatea  between  0  ana  0.5  ope  fro.  the  solutron  at 
rr^  vet  we  expeot  tne  correct  aofotione  to  beco.e  .ore 

the  array  aperature  .easured  in  wavelengths  heco.es 
erratic  as  ar  y  P  frequencies, 

11  Tf  the  optimum  filt'^rs  wer^^ 

smaller.  If  tne  p  absent,  as  seems  likely  from 

and  if  the  second  higher  mode  were  ah  the  5  x  3  and 

.  ™  „n  Figure  7:  then  the  outputs  from  tne 

the  bottom  trace  on  Figure  /, 

the  5  X  4  filters  would  match. 
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IV.  CONCLUSIONS 


Qualitative  experierce  has  disclosed  the  following  features 
of  our  optimum  filter  solutions; 

1.  These  filters  arc  indeed  orthogonal  and  separate  the 
noise  md  signal  modes  as  planned. 

2.  These  optimum  filters  extend  for  400  points  in  time 
since  the  frequency  interval  chosen  was  .05  cps. 

3.  The  gain  on  the  synthetic  examples  was  greater  than 
20  db  with  solutions  restricted  to  0.5  to  2.0  frequency  range. 

More  resolution  in  frequency  could  increase  this  figure. 

4.  The  optimum  filters  are  zero  phase  shift  filters. 
Therefore,  they  are  algebraically  additive. 

5.  The  optimum  filters  become  larger  in  gains  (both  positive 
and  negative)  as  degrees  of  freedom  decrease. 

6.  The  optimum  filters  become  larger  in  gains  (both  positive 
ard  negative)  as  the  aperature  of  array  goes  down.  Thus,  for  low 
frequencies  the  solutions  tend  to  become  unstable. 

7.  Extra  modes  not  considered  and  errors  in  'assumptions 
cause  errors  in  the  output.  Errors  in  assumptions  can  include 

an  incorrect  well  log,  non-parallel  layering  in  surrounding  medium, 
and  incorrect  calibrations  of  seismometers. 

8.  Extra  degrees  of  freedom  are  needed  to  cut  gains  of 
optimum  filters  and  make  optimum  solutions  more  tolerant  of  errors. 

9.  Extra  degrees  of  freedom  are  best  obtained  by  increasing 
the  number  of  seismometers  in  the  vertical  array- 

10-  More  stable  solutions  (i.e.,  optimum  filters  with  lower 
gain)  will  be  obtained  from  the  deeper  vertical  arrays  which  have 
the  seismov/.ieters  distributed  rather  uniformly  throughout  the  array. 


we  will  evaluate  these  properties  of  vertical  =»  rays 
quantitatively  as  the  -tudy  progresses,  we  will  ar  .o  add 

deghosting  to  the  optimum  filter  programs. 

The  next  vertical  array  of  Interest  is  the  one  in  the 
Grapevine  well,  GVTX.  This  array  has  six  seismometers  down¬ 
hole  plus  a  surface  instrument  all  recorded  on  the  same  tape, 
in  addition  it  is  located  in  a  sedimentary  basi..  rather  than 
the  Rocky  Mountains  so  tie  assumption  of  plane,  horisontal 
layering  may  be  better  satisfied. 
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